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We study the effects of symmetry breaking on non-Markovian dynamics in various system-bath 
arrangements. It is shown that by breaking certain symmetries features signaling non-Markovian 
time evolution disappear within a finite time t g . We demonstrate numerically that the scaling of 
t g with the symmetry breaking strength is different for various types of symmetry. We provide 
a mathematical explanation for these differences related to the spectrum of the total system-bath 
Hamiltonian and provide arguments that the scaling properties of t g should be universal. 



Introduction 

Non-Markovian dynamics occurs in systems connected 
to a bath when information travels from the bath into 
the system. The name non-Markovian arises from the 
fact that these systems cannot be described by Marko- 
vian master equations [TJ [SJ. A system which is non- 
trivially joined to a finite bath will always be necessarily 
non-Markovian because any information which leaves the 
system must return [3]. In a generic finite system, the 
recurrence time is expected to scale exponentially with 
the size of the Hilbert space(4]. For many systems with 
a moderately sized finite bath, it may therefore be the 
case that significant non-Markovian events are only ob- 
servable on timescales much longer than those of other 
relevant dynamics. 

In order to define and quantify non-Markovianity in 
quantum systems different methods have been devel- 
oped. For example, one can examine a distance mea- 
sure between the actual evolution and the best Marko- 
vian model of a system bath arrangement. However, this 
approach is rarely used because of the difficulty of find- 
ing an optimum Markovian approximation |5j. There are 
other methods which involve examining the mathemat- 
ical structure of the map defining the time evolution of 
the system [BHH]- Other techniques which have been pro- 
posed are based on determining the minimum amount 
of noise required to make the dynamics Markovian [5] or 
studying the evolution of entanglement for a system pre- 
pared initially in a maximally entangled stateJS]. 

Here we choose a method of studying the degree of 
non-Markovianity similar to the one proposed in and 
further examined in [T2l [T5] . Its central idea is to exam- 
ine the trace distance between two different initial sys- 
tem states under time evolution (for examples of other 
methods and techniques, both analytic and numerical, 
see [21 15 HTU1 IT5H^T| V The trace distance is defined as 

= ^Tr(yW, \M) - Ps(t, \%W, (l) 

where |^o) and \i/j' ) are two different initial states 
which are product states of the system and the bath with 
the same bath states but orthogonal system states. The 
system density matrix at a time t is defined as 



p s (t,\ij )) = Tx bath (m))^(t)\), (2) 
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where H is the overall system-bath Hamiltonian. A 
system is necessarily behaving in a non-Markovian fash- 
ion if the trace distance increases in time because of its 
interpretation as a measure of distinguishability. In gen- 
eral, the maximum of (t) over all possible initial 
states has to be calculated (see Ref. |TTj). However, in 
the case where the system has only a few degrees of free- 
dom, a single pair of orthogonal states can still provide 
valuable insights into the overall non-Markovian behavior 
of the system. Furthermore, examining non-Markovian 
behavior of a system by studying only two orthogonal 
initial states has the advantage that in numerical studies 
it is often impractical to maximize over all possible initial 
state pairs. 

Non-Markovian systems have been studied intensively 
both theoretically and experimentally in recent years. On 
the experimental side the transition between Markovian 
and non-Markovian dynamics has been probed in opti- 
cal setups [2"5l \M\ . These experiments have been able to 
go from one regime to the other in a controllable way, 
thus allowing them to control the flow of information be- 
tween system and bath. Non-Markovian effects have also 
been observed in solid state systems, for example with 
electron transport though a quantum dot |25| . Theo- 
retically, recent research in the field of non-Markovian 
dynamics has been devoted to e.g. the study of memory 
effects^ [201 IS]> including examining the effect of bath 
memory on the equilibrium state of a system after evolv- 
ing for a long timelJD], the influence of a chaotic environ- 
ment on non-Markovian dynamics |22| or the connection 
between non-Markovianty and entanglement [26J. Still, 
our knowledge on the question how the internal struc- 
ture of the bath influences the dynamics of the system is 
incomplete. 

In this paper we demonstrate that non-Markovian dy- 
namics can be controlled by breaking symmetries (e.g. 
parity) of the complete system-bath Hamiltonian. We 
show that non-Markovian events in the temporal evolu- 
tion of the trace distance disappear on a certain time 
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scale which depends on the strength of the symmetry 
breaking and give spectral arguments to show that the 
scaling behavior we observe is universal. 

This paper is structured as follows. In Sec. [T]we explain 
the two system bath Hamiltonians studied in this paper, 
as well as the initial states and why they are chosen. In 
Sec. [II] we look at the dynamics of a tight binding torus, 
which has combined system-bath symmetry. In Sec. |III| 
we look at a 2 level qubit system coupled to a completely 
connected tight binding bath. In this case we break a 
symmetry which exists only in the bath. This is followed 
by a brief summary of conclusions. More detailed calcu- 
lations are given in Sec. |IV| 



I. SETUP 

The first setup we study is a two-dimensional tight 
binding lattice whose dynamics is governed by the Hamil- 
tonian 
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Figure 1: Illustration of the two models studied in this pa- 
per a) system: 2 tight binding sites, bath: surrounding 2 
dimensional tight-binding lattice b) system: single spin-| in 
a magnetic field, bath: fully connected tight binding graph. 
Quench: system-bath coupling turned on at t=0. 
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where Ci(c]) is the annihilation (creation) operator on 
site i and N is the number of sites in the lattice. The 
< i,j > notation indicates that the sum is performed 
only over adjacent sites, ry is a symmetric matrix, i.e. 
r ij = r ji w ifh f ne additional property |ry| < 1, whose 
elements are chosen in a uniform random fashion, g is 
a positive real number measuring the "strength" of the 
symmetry breaking. In the case of g = we recover the 
Hamiltonian of a particle which can move freely between 
the sites of the lattice. For finite values of g the matrix 
elements of the Hamiltonian are altered randomly, which 
corresponds to altering the hopping rates between the 
sites. We initialize the system in the ground state of 
the Hamiltonian where 2 sites are disconnected. In the 
following these sites are referred to as the "system" and 
the remaining sites are referred to as the "bath". At t = 
the system and the bath are joined in a "quench" process. 
Fig. [T^l) shows the lattice and the quench. Evidently, 
the two eigenstates of the two site system Hamiltonian 
without any connection to the bath are a good choice of 
initial system states for the trace distance analysis. 

The second setup we consider is a two- level (spin |) 
system coupled to a bath which is a fully connected 
graph. We have chosen a fully connected graph as the 
bath Hamiltonian because it has a very high degree of 
symmetry (symmetry under any permutation of sites) 
and therefore is a natural choice for examining the ef- 
fects of symmetry breaking. The corresponding complete 
system-bath Hamiltonian is given by 
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which is of the canonical form H = _ff. s 
Hh. 



sys. yy '-bath ' 

sys. *x "bath t- k H^p ® HlZl- As before, r' pq and r %] 
are real symmetric matrices, (r' = r' . rjj = Tji) with 



H conn (g,k,N,m,r, r') 



'-bath- 



(4) 



' pq ~ r qp' ri '-> ~ ' i % > 
the property \r' pq \ < 1 , \nj\ < 1. The elements of r' pq 

and Tij are selected in a uniform random way. k mea- 
sures the strength of the coupling between the system 
and the bath. For k = the system and the bath are 
completely decoupled. As k is increased a coupling is 
generated between the direction of the spin and the tun- 
nelling constants in the subset of sites 1 through m of 
the bath. We consider the case where < k < 1. De- 
pending on the sign of r^this coupling leads to either 
an increase or decrease in the hopping strength on each 
bond in the subset. If the direction of the system spin 
is switched, the hopping probabilities which were more 
favorable will become less favorable and vice versa. As 
in the first case g is a real number which measures the 
strength of the symmetry breaking. Fig. [T]d) shows a 
sketch of this system-bath arrangement. For this Hamil- 
tonian the spin is prepared in the +z (-z) direction and 
the bath is initiated with the lower energy quarter (25%) 
of its k — eigenstates filled with non-interacting par- 
ticles. By means of this choice we can maximize the 
amount of information exchange between the system and 
the bath (for details see Sec. |IV B ). At t = a quench is 
performed when the spin-bath coupling k is turned on. 
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Figure 2: a) Trace distance versus time for two 10x10 lat- 
tices, (red) circles are for a torus and (cyan) diamonds are 
for a strip, b) Trace distance versus time for a 10x10 torus at 
half filling for various strengths of randomness g. g=0 (dotted 
blue line), g=0.02 (dashed red line), g=0.05 (solid green line), 
g=0.1 (dot dashed purple line), c) Timescale at which recur- 
rences disappear, t g versus g for 10x10 torus at half filling, 
upper solid noisy (blue) line is the mean of i 9 and the dashed 
(green) line is a linear fit. The lower solid (red) line is the 
median of t g and the dot-dashed (black) line the correspond- 
ing linear fit. d) Non-Markovianity measure ,yV(t, At; D) with 
At — 1000 (blue) Xes are the mean of the measure for 400 re- 
alizations from t = 0, (red) squares are the measure performed 
on the mean trace distance at times from t = 0, (green) as- 
terisks are the mean of the measure for times from t — 10 s 
and (black) diamonds are the measure of the mean at times 
from t = 10 8 . All data in this plot are for small randomness, 
5 = 10" 6 . 



II. TWO DIMENSIONAL TIGHT-BINDING 
LATTICE WITH COMBINED SYSTEM-BATH 
SYMMETRY 



Let us now discuss how the trace distance (t) 
evolves in time for the tight binding lattice. For a dis- 
cussion concerning the numerical details, such as how we 
calculate the density matrix, we refer to the supplemen- 
IV Al Fig. [2ji) shows D^ {t) for 



tal material (Sec. 

two different boundary conditions applied to the Hamil- 
tonian H2D given by Eq. [3] For periodic boundary con- 
ditions in x (x and y) we get a strip (torus). Since we 
monitor the temporal evolution of two orthogonal ini- 
tial states the trace distance is one in the beginning of 
the dynamics when system and bath are joined together, 
that is D^>(t = 0) = 1. Afterwards D^>(t) decreases 
rapidly in a short time, which is not visible in Fig. [2J1,) 
and then continues to fluctuate around a small value 
D^ifi'it — 0) « 0.35. In the case of the torus, which has 



the greater amount of symmetry of the two geometries, 
these fluctuations of the trace distance D^> (t) are inter- 
rupted by pronounced peaks. By performing a thorough 
analysis of the dynamics of the setup at these points in 
time we have found that these peaks can be explained by 
partial reconstructions of states which are similar to the 
initial states having a trace distance of one. In the follow- 
ing we elucidate this mechanism of state reconstruction. 

For the torus the outgoing wavefunction from the ini- 
tial state which is localized on the system sites propa- 
gates freely in both directions when the connections be- 
tween the system and the bath are turned on. Eventually, 
it returns coherently to its initial position, reconstruct- 
ing a wavefunction which is mostly isolated within the 
system. Similarly, the wavefunctions of the other states 
which initially have zero amplitude in the system recon- 
struct states with an approximately similar shape. On 
the contrary, for the strip geometry the wavefunctions 
are scattered off the open boundary such that they do not 
reconstruct coherently when they return. Formally, the 
torus has s.'Ln symmetry group whereas the strip 

only has the cyclic symmetry group Zjyj35], The larger 
symmetry group of the torus allows the wavefunction to 
preserve more of its shape while propagating. Therefore 
the torus allows for reconstructions which do not occur 
for the strip. Note that these reconstructions are not the 
same as recurrences because any finite closed quantum 
system must eventually return to a state arbitrarily close 
to its initial state [3] • Recurrences occur in the strip as 
well as any other geometry we could consider, but they 
occur at a much longer timescale than the partial recon- 
structions. 

The coherence of the returning waves which partially 
reconstruct the initial state can be destroyed over time 
if randomness is introduced to the Hamiltonian. Fig. 
2b) shows the temporal evolution of the trace distance 



D^'{t) in the case of the torus-geometry for different 
values of g, the strength of the random term. As the mag- 
nification in Fig. |2|3) reveals the height of reconstruction 
peaks decreases with increasing g. In order to analyze 
this behavior we have studied how the time it takes for 
these peaks to disappear t g depends on the strength of 
the randomness g. Fig. pp) shows t g as a function of g in 
a double-logarithmic plot. As we can see from Fig. [2J1) 
t g depends on g in an inverse linear fashion, i.e. 



t g <x 



(•5) 



In the following we derive an intuitive understanding 
for this behavior. Evidently, the Hamiltonian of the com- 
plete system-bath arrangement has a discrete transla- 
tional symmetry which leads to level-degeneracy in the 
energy spectrum. The addition of randomness breaks 
this symmetry and therefore leads to level splitting. For 
small randomness the gaps, which are created by the level 
splitting, are linearly proportional to g. The timescale 
associated with these gaps is simply their inverse which 
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yields the dependency of t g on g given by Eq. (|5j). The 
disappearance of the reconstruction by introducing ran- 
domness to the system-bath Hamiltonian peaks can be 
thought of as an equilibration process. We find that away 
from the reconstruction peaks the system density matrix 
is close to the fully equilibrated state for both initial con- 
ditions. The reconstructions drive the system away from 
this state, and therefore the system becomes more equili- 
brated as they are destroyed by symmetry breaking. The 
timescale t g with which these peaks disappear can there- 
fore be thought of as a timescale for the equilibration of 
the system. (We also note that in Fig. ^p) the mean of 
t g is quite noisy and is much greater than the median. 
This indicates the presence of rare events where t g is very 
large.) 

Finally, we examine the effect of the reconstructions 
on the non-Markovianity of the system. To this end we 
calculate the total increase in trace distance over a given 
time period defined as 

t+At 

jr{tM,D)= J dr[ ^l)e(^M)], (6) 

where is the Heavyside function and D^i (t) is de- 
fined in Eq. [I] jY{t, At, D) is a measure of the informa- 
tion flowing from the bath into the system which allows 
us to distinguish between tp and ip'. It can therefore 
be thought of as the information allowing us to distin- 
guish between ip and ip' which would necessarily be lost 
by a Markovian approximation of the bath. Any time 
distinguishbility between two arbitrary states is increas- 
ing the behavior of the system can be said to be non- 
MarkovianfTT]. Fig. [2ji) shows jV(t, At, D) as a function 
of the number of system sites N S i tes at different points 
in time t with a fixed At. At early times the measure is 
the same for all realizations because the symmetry is not 
broken yet so the mean of the measure is the same as the 
measure of the mean. We see that with increasing time 
jV{t, At, D) decreases systematically for all system sizes. 
Consequently, the evolution is less non-Markovian after 
the reconstruction peaks have disappeared. Furthermore, 
the average over many runs shows asymptotically almost 
no non-Markovianity, indicating that after the symmetry 
is broken there are no non-Markovian features which are 
independent of the random perturbations to the lattice. 
Since the coupling between the system and bath is quite 
strong, there is no reason to expect Markovianity even 
asymptotically and in the large bath limit . 

Moreover, Fig. |2jl) reveals that as the system 
size N s ites becomes larger the early and late non- 
Markovianity measures approach each other. An expla- 
nation for this is that as the system is made larger there 
is more time between the reconstructions. Thus the ef- 
fect of the reconstructions is hidden by small fluctuations 
which transport information in and out of the system at 
all times. 
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Figure 3: a) Trace distance between two initially orthogonal 
states versus time for a bath consisting of 25 particles on 100 
sites (N=100) with k=l and m=20. Pluses (green) correspond 
to g = 0, circles (red) to g — 10 -3 , and diamonds (black) 
to g = 10~\ b) t g versus the strength of the randomness g 
respectively the coupling parameter between system and bath 
k. The solid (blue) line is t g versus g for the same parameters 
as (a), the dot-dashed (magenta) line is the same but with 
k — 10~ 3 . The dashed (red) line is t g versus k for g=l and 
all other parameters the same as (a), the dotted (green) line 
is the same but with g — 10~ 3 . c) D^>(t) for various system 
sizes with k = 1, m = ^ , and quarter filling. Solid (blue) line 
is for TV = 100, dashed (green) TV = 200, dot-dashed (red) 
TV = 300, dotted (cyan) TV = 400, dashed (purple) TV = 500 
and solid (gold) TV = 600. All plots are averaged over 400 
realizations of the system-bath arrangement, d) Standard 
deviation of mean trace distance averaged over 400 samples 
for g = 10~ 4 . Standard deviation of mean (red asterisks) 
and mean of standard deviation (green Xes) for times from 
t = 10 5 to t = 10 7 , standard deviation of mean (black circles) 
and mean of standard deviation (blue squares), for times from 
t = 10 9 to t = 10 11 . 



III. FULLY CONNECTED GRAPH COUPLED 
TO 2 LEVEL SYSTEM 

Let us now turn to the second setup which is a qubit 
coupled with a fully connected tight binding graph. Fig. 
|3^i) shows the temporal evolution of the trace distance 
D^j/j'(t) for different values of the randomness g. For 
g = we observe fluctuations of D^r(t) around 0.8 after 
a rapid decrease from the initial value D^>(t = 0) = 1. 
When the randomness is large, g — 10 _1 , the fluctuations 
occur at a much smaller value, that is D^i(t) ~ 0.2. 
For moderate randomness, g = 10~ 3 , the trace distance 
slowly decreases in time until it continues to fluctuate 
around an intermediate value between 0.2 and 0.8. Ac- 
cordingly, it can be said that the system equilibrates to 
states with a much larger trace distance for g = com- 
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pared to the case when g is finite. Since the trace distance 
tells us about distinguishability and the information con- 
tent of the system, it is considered as a measure of equi- 
libration. In other words one can say that the unbroken 
permutation symmetry for g = prevents full equilibra- 
tion for which we give a physical explanation in the next 
paragraph. 

For g — the high degree of symmetry leads to many 
completely localized eigenstates in the part of the bath 
which is not connected to the system. Consequently, 
these parts of the wavefunction cannot contribute to the 
equilibration process. Hence for g = the system will 
never fully equilibrate. When a small randomness (g ^ 0) 
is added, these eigenstates acquire a finite amplitude to 
be in the part of the bath connected to the system and 
thus the system will eventually equilibrate, (see Sec. 
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Contrary to the first setup we observe in Fig. |3[ 
no peaks or systematic increases of the trace distance 
D^^i (t) which would correspond to a flow of information 
from the bath into the system. Consequently, one could 
draw the conclusion that the second setup shows no sys- 
tematic non-Markovian dynamics which do not depend 
on the details of r' pq and r^. However, this is not true 
and we will return to this problem at a later point in the 
manuscript . 

Let us now ask how the strength of the random term 
g affects the equilibration time t g which we define as the 
time when the trace distance drops below 0.4. In Fig. 
3j}) t g versus g is shown in a double logarithmic plot . As 
we can see the equilibration time scales as 



1 

t g OC — . 

g 2 



(7) 



as long as g is sufficiently small compared to the cou- 
pling parameter k. Furthermore, we observe that if g is 
large compared to k the equilibration time t g does not 
depend on g. 

The physical reason for the different dependencies of 
the equilibration time t g on the strength of the random- 
ness g for the two studied setups (compare Eqs. [5] and 
[7]) can be traced back to the spectrum of the system 
bath Hamiltonian. Both Hamiltonians show a level split- 
ting which is proportional to g. Nevertheless, this level 
splitting leaves a gap degeneracy in the case of the fully 
connected bath which is effectively the same as a level de- 
generacy once we trace out the bath degrees of freedom. 
A demonstration that this happens whenever there is a 
symmetry in the bath which is not broken by coupling to 
the system can be found in Sec. |IV C| 

As the next step we keep the randomness fixed and 
study how the equilibration time t g depends on the 
strength of the coupling between the system and the bath 
k. Fig. [3]b) shows that t g scales exactly as before in an 
inverse quadratic way, i.e. 
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which we are going to explain in the following. In this 
case the tensor product nature of the k=0 Hamiltonian 
produces a gap degeneracy. The tensor product structure 
of the coupling term in Eq. [4] means that this gap de- 
generacy cannot be broken by direct level splitting, but 
instead must be broken by level repulsion, for more de- 
tails see Sec. ITVDl 

Now we return to question whether the dynamics 
of the second setup possesses hallmarks of systematic 
non-Markovianty or not, which has come up since Fig. 
[3^,) shows no sytematic increase of the trace distance 
D^>(t). Fig. [3k) shows an average of D^^i(t) over 400 
different realizations of r' pq versus the bath size whereas 
the fraction S is fixed and we choose k — 1, g — 0. Evi- 
dently, we observe systematic non-Markovian oscillations 
in Fig. p3p) which are masked by random fluctuations in 
Fig. [3k^ Indeed these oscillations are systematic in the 
sense that they are not destroyed by averaging over dif- 
ferent realizations. Increasing the bath size appears to 
actually increase the strength of the systematic oscilla- 
tions. These oscillations are caused by internal bath dy- 
namics where the bath sites which are not connected to 
the system effectively act as a single site, detailed calcu- 
lations of how this happens can be found in Sec. |IV E| 

After the system has fully equilibrated the non- 
Markovian oscillations of the trace distance disappear as 
we can see from Fig. [3ji) in which we examine the stan- 
dard deviation of the trace distance at early (after the ini- 
tial transient but before the symmetry is broken) and late 
(after the symmetry is broken) times. This can be seen 
because at late times the standard deviation of the mean 
over different choices of r\j and r' pq drops to nearly zero. 
The drop to nearly zero means that the systematic os- 
cillations have disappeared because the mean value over 
different realizations hardly fluctuates at all. At these 
late times we can see that there are fairly strong fluctu- 
ations in D^,i (t) , these fluctuations disappear however 
when we average over different realizations. The ana- 
logue of Fig. [3]:) at late rather than early times would 
look like a flat line and not contain the sinusoidal oscil- 
lations we see at early times. 

By examining the standard deviation of the mean trace 
distance compared to the mean of the standard deviation 
we can determine how sensitive the non-Markovian be- 
havior is to the specific values of the randomly selected 
term r' (and for late times also r^). As Fig. [3tl) shows 
the standard deviation of the mean of trace distance at 
early times increases and is approached by the mean of 
the standard deviation for large system size. The close- 
ness of these two numbers indicate that for large bath 
sizes the trace distance at early times is dominated by 
behaviors which are the same for most choices of r' pq and 
do not depend much on the details of the random modu- 
lation of bond strengths. In contrast the standard devia- 
tion of the mean trace distance for late times is much less 
than the mean of the standard deviation for all system 
sizes, indicating that asymptotically the trace distance 
is dominated by random fluctuations strongly dependent 



6 



on the details of nj and r' pq . 

In Fig. |3ji) we examine the mean of the standard devia- 
tion of the trace distance over different random selections 
of Tij and r' pq , we also examine the standard deviation of 
the mean value of trace distance. The mean of the stan- 
dard deviation tells us about the average variability of 
the trace distance before and after the symmetry is bro- 
ken. This mean decreases with system size for late times 
and eventually becomes less than the value for early times 
for a bath with 800 total sites. This decrease with size is 
caused by a decrease in non-Markovian fluctuations asso- 
ciated with finite size effects as the bath is made bigger. 
In contrast the mean of the standard deviation for early 
times increases with system size as the trace distance os- 
cillations shown in Fig. [3t) get stronger. 



IV. SUPPLEMENTAL MATERIAL 

A. Details of density matrix construction for 2 
dimensional graph 

In this case we are considering a bosonic system and 
producing a reduced density matrix which combines the 
matrix elements for measuring a single particle on a site 
with those of measuring any finite (non-zero) number of 
particles on that site. For a two site system the density 
matrix for an individual wavefunction can be written as 



WM\i>) (9) 

where 



Conclusions 



In this work it has been shown that non-Markovian dy- 
namics can be evoked by symmetry breaking. We have 
examined two examples of atypical equilibration which 
involve highly non-Markovian behavior for a period of 
time. In both of these examples, internal bath degrees 
of freedom play an important role. We have shown that 
these non-Markovian features are related to symmetries 
in the Hamiltonian, and that breaking these symmetries 
determines the timescales of the disappearance. Fur- 
thermore we have demonstrated that the effect of the 
symmetry breaking process is fundamentally different for 
symmetries of the overall system-bath Hamiltonian. com- 
pared to the breaking of symmetries which exist only in 
the bath, but are not broken by coupling to the system. 
There is a difference in the scaling of the equilibration 
timescale with the strength of the symmetry breaking. 
This difference arises from the fact that for the break- 
ing of a combined system-bath symmetry the relevant 
equilibration timescale is related to the breaking of level 
degeneracy, whereas in the case where a bath symme- 
try is broken the equilibration timescale relates to the 
breaking of gap degeneracies. Because of the underlying 
spectral cause, these types of scaling are expected to be 
universal. 
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where Si.j is a Kronecker delta and <0,0> denotes the 
set of states where neither site is occupied, <1,0> denotes 
the state where site 1 is occupied and <0,1> denotes 
the state where site 2 is occupied. These single particle 
density matrices can then be combined into a 4x4 total 
density matrix by taking the product of density matrix 
elements to produce a combined density matrix using the 
following rule 
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where B is a 4x4x4 tensor defined as 
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These rules can be applied repeatedly to build up a full 
multi-particle density matrix. In this density matrix the 
fourth diagonal element corresponds to no particles on 
either site, the third diagonal element corresponds to a 
non-zero number of particles on site 2, but zero particles 
on site 1, the second corresponds to none on site 2 but a 
non-zero number on site 1, and the first corresponds to 
a non-zero number on both sites. 
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e 

Figure 4: Estimate of the long time mean of trace distance 
versus 9, the angle between ipo and the x-axis with ip' Q chosen 
as the orthogonal state to tf> . 



B. Estimate for long time mean of trace distance 

We now argue how the relevant non-Markovian dynam- 
ics of this system can be illuminated by observing only a 
single pair of orthogonal states. For the case of the spin 
system the choice of the two orthogonal initial states is 
not as obvious as it was for the torus. As in the case of 2d 
geometries, every possible choice of initial state in of the 
system has a single uniquej^S] orthogonal partner. How- 
ever, the best choice of initial state is not clear a priori. 
In the following we explain the choice of initial states for 
this setup. H sys rotates spin vectors around the x-axis. 
Therefore, a rotation of the initial pair of states about 
this axis amounts to connecting the system to the bath 
at a different time. Fig. ^demonstrates that the average 
trace distance should be smallest when ip is in the y-z 
plane. We can understand this intuitively because ini- 
tial states pointing in the ±x direction are eigenstates of 
H sys . ® hath + hys. ® H bath , whereas states in the y-z 
plane are only eigenstates of l sys . (8 Hbath- We want to 
examine states in which the maximum information trav- 
els into and out of the system, so we choose states in the 
y-z plane. 

The estimate of the long time mean is calculated 
by totally dephasing (i.e. setting all off diagonal el- 
ements in the energy basis to zero) the total system- 
bath density matrix for the states i/jq and i/j' and cal- 
culating the trace distance between the resulting sys- 
tem density matrices. Let us call the totally de- 
phased matrix p^, for any operator A, Trf.ApooGV'o))] = 
f™ dTTr[Aexp(-iH r)|V>oX^o|{exp(iEf r))], the long 
time mean of A. The trace distance can be though of as 
the maximum possible distinguishibility between systems 
with the two initial states using only measurements in the 
system, i.e. D^(t) = maxA(Tr[ATr bat h[p(t, |^o})]]) - 
Tr[ATr oat h[p{t,\il)' ))})) where the eigenvalues of A are 
constrained to lie between 1 and 0. We therefore 
have for the trace distance between Tibathlfioo {l^o})] & n d 



^bath[Poo{\lp' ))}, 

max{Tr[ATr bath [pMo))]}) - Ti[ATrbath[pMo))}]) 
= max(Tr[^4(Tr ha i, i [p oo (|-0 o ))] - Tr ba th[Poo(\4'o))})} 
= max(Tr[^l(Tr ha4 , l [p 00 (|i/'o>) - Poo{\i>' ))]] 

1 f* 

= max( lim - / dr(Tr[A(exp(— iH t)|^oaW exp(i7J r)) 
A t^oo t J T=0 

- exp(-ifrr)|^M|exp(»^r))])) 

i r* 

< lim - / dr(max(Tr[A(exp(-ii? T)\ip Q )(ip \esq>(iH t)) 

~ t^oo t J T=0 A 

- exp(-iffr)|^XVo|exp(iffr))])) 

1 f* 

= lim - / iJtDjjmM. 

t^oo t J T=0 

The trace distance between Trf, a t? i [poo(|V , o))] and 
r ^bath[Poo(\' t Po})] gives the maximum average distin- 
guishability between p s (t, \ipo}) an d p s (t, \ip'o)) assuming 
measurements are performed in the same basis at all 
times. The trace distance between the totally dephased 
states provides a lower bound for, and estimate of the 
average trace distance between p s (t, \ipo)) and p s (t, \ip' }) 
which is the maximum average distinguishability if mea- 
surements can be made in a different basis at every time. 

C. Scaling arguments for t g versus g for bath 
symmetry being broken 

An explanation for why t g scales like p- rather than ^ 
can be obtained by examining the energy spectrum of the 
Hamiltonian at different values of g. As we can see from 
Fig. |5^l) , the spectrum for g=0 is highly degenerate. As 
g is increased, the degeneracy is broken, however there 
will be pairs of states within the two manifolds which 
change energy with the same slope, as we can see from 
Fig. |5|l). The gap between the energies of the states in 
these pairs will therefore remain the same until the slope 
of the energy is changed by an avoided crossing. The 
relevant process for equilibrating the system is not the 
splitting of the level degeneracy but rather the splitting 
of a gap degeneracy caused by avoided crossings. Because 
the slope is the same for each pair, the gap must change in 
a manner proportional to g 2 rather than g. We can verify 
that the gap increases proportionally to g 2 by looking at 
the double logarithmic plot of gap versus g shown in Fig. 

^' . . r. 

The two degenerate manifolds in ph) correspond to 

states where the particle in the bath is isolated in the 
sites which are not connected to the system. These two 
manifolds correspond to the two eigenstates of the sys- 
tem Hamiltonian. The total Hamiltonian can therefore 
be written in the form 

sym) © H nosy m 

(ii) 
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Figure 5: a) Energy spectrum of the Hamiltonian for vari- 
ous strengths of g, dot-dashed (yellow) lines are eigenstates 
which are not initially part of a degenerate manifold of states 
for g=0, dotted (red) lines are states which are part of a de- 
generate manifold for g=0. A pair of eigenstates which are 
the same in the bath for g=0 are shown as sold (blue) lines. 
The slope of the energy of these two states is the same for 
small g but is changed by avoided crossings for larger g. b) 
Double logarithmic plot of the difference from the initial gap 
between pairs of states within the degenerate manifold which 
are the same in the bath for g=0 versus g. 



gap between states within the different degenerate man- 
ifolds which is determined uniquely by n and m. We 
can now trace out the bath, yielding the system density 
matrix, 



Pnm{t) = y^(exp(-zA nm )( A qqlkmn ) (15) 

q l,k£sym 

+ 2J A qq i kmn exp(-iAi k t)) 

I or k(£sym 



Notice that all of the terms with Ay ^ do not con- 
tribute to the final reduced density matrix, and therefore 
the sector of the bath containing H sp ut cannot dephase 
the system. If X m ix 7^ than level repulsion can occur 
with energy levels outside of the subspace as shown in 
Fig. [5] these level repulsions break the gap degeneracy 
and therefore destroy the tensor product structure in Eq. 
[l2| As we can see from Fig. ^p) the gaps are widened pro- 
portionally to g 2 rather than g, therefore the timescale for 
the broken symmetry to equilibrate the system is propor- 
tional to ■ The arguments given here will work for any 
case where a bath symmetry is broken, therefore for small 
g, t g oc for a generic bath symmetry being broken [37], 



This is the generic form of any system bath arrange- 
ment with a bath symmetry. In general we can always 



write H, 



break 



"split {H S plit 



H 2 ) + X mix H mix where 



H S piit spans the same subspace as l sym in Eq jll| The 
Hamiltonian is now of the form 



H — (H' sys (g) l sym + g X split Isys <£> H split) (12) 
®i,H noS ym "F Q X S plit H2) "F Q Xmix lsys ® H m i x . 

We can now construct a basis in the system such that 
In) (n|, if'„J = V|n) and a basis in the bath such that 



I*) = 



1 \i) £ deg. 



where deg. refers to the 



(J otherwise 

degenerate subspace created by the bath symmetry. For 
X m ix — the density matrix of a system with an arbitrary 
initial state l^o) can then be written as 

pijmn(t) = exp(— iA mn t) (13) 
• ( AijUemneaqpi-iAijt)) 

+ ^ Aiji kmn exp(-iAi k t), (14) 

I or k^sym 

{{m\®{i\)\l) (i> \k) (fc|(|j>®|n)), 
Ei-E k 



^■ijlkmn 



Where \l) and |fc) are eigenstates of the overall Hamil- 
tonian and A mn = (m \ H' sys \ m) — (n \ H' sys \ n) is the 



D. Scaling arguments for t g versus k 



We ca n ma ke very similar arguments to those given 
in Sec. IV C for why t 9 oc p. First let us consider 
the Hamiltonian of a system coupled to a bath in the 
canonical form 



H = H 



sys. 



Lbath- 



<-sys. 



bath ' 



TjCOUp 

> n bath ■ 



(16) 



Note that Hamiltonian given in Eq. [4] is of this form. 
We now note that we can decompose H^ p (g> H^°^ — 

where [D c ™P,H sys .] = 



recoup 
sys 



DbaZ+R C s°y7< 



rycoup 
U bath 



and [Dba^th , Hbath] — 0. We now note that we can 



write H s 

TDCOUp 

sys. 



i-bath + If 



H, 



bath 



r>coup\ 
^bath > 



HDTy 



bath 



R Zth where R ii k ) 



cannot equilibrate the bath by itself because [D^ p ® 

DZLHsys. ® hath] = and [D™ P ® DZlAsys.® 

Hbath] = 0. In other words, a gap degneracy will persist 
in Hi (k) for any value of k. The gap degeneracy is broken 
by the term k Rf^ p ® B%%£, however <g> R^* can 

always be chosen such that it contains only off diagonal 
elements in the energy basis of H\(k = 0). The off di- 
agonal elements can only change the spectrum though 
level repulsion rather than direct level splitting. The 
deviation from gap degeneracy will therefore scale like 
fc 2 rather than k, and as a result we observe t g oc rj- 
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E. Eigenstates and eigenvalues of fully connected 
graph 

Consider the Hamiltonian of a fully connected tight 
binding graph with equal coupling between all sites and 
zero on-site potential, H c = X)fj c l c j- This Hamiltonian 
has a N-l fold degenerate ground state manifold that 
spans all states for which the total amplitude sums to 
zero. i.e. ipi = 0. This can be shown in the case of zero 
on-site potential by observing that (H c ip)i = J^j^i^j- 
Therefore in the case where Yli^i = 0) we see t na * 
(-HcVOi = Ylj V'j — — —ipi) these states therefore are 
eigenstates with eigenvalue -1. In the case where ipi = 

*i = CJ vf we have ( H ^)i = = (N- 

this state is therefore an eigenstate with an eigenvalue of 

N-l. 

Now let us now consider a Hamiltonian of the 
form H ' comb . = H c (B H a + H join , where H Join = 
53^1 1 Sj3v+i c l c j equally couples all N sites in H c to 



another group of sites, H a is an arbitrary Hamiltonian 
of size M. We now note that for an initial state ip tra P 
such that ^T ap = and = V i > iV. We 

can show that ip tra P is an eigenstate of H c © H a be- 
cause it has no support for i > N, and the part for 
i < N is an eigenstate of H c with an eigenvalue of 
(N-l). We now also see that (Hj oin ->jj trap ) i < :N = and 

{H join ip trap ) i>N = J2iLi ^T aV = °! therefore H comb will 
have an N-l fold degenerate manifold of states which only 
have finite amplitudes in the first N sites and will have 
an energy of N-l. Wavefunctions in these states are ef- 
fectively trapped in the first N sites. 

We know by orthogonality of eigenstates, that all other 
eigenstates of H C omb must have tpi = ipj = ^ cc for all i 
and j less then N+l. This constraint means that the in- 
ternal degrees of freedom of the first N sites do not play a 
role in the eigenstates outside of the manifold mentioned 
in the previous paragraph. The first N sites can therefore 
be modeled as a single site. 
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